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Abstract 

The vacuum polarization functions n(g 2 ) of charged and neutral gauge bosons which 
arise from top and bottom quark loops lead to important shifts in relations between 
electroweak parameters which can be measured with ever-increasing precision. The 
large mass of the top quark allows approximation of these functions through the 
first two terms of an expansion in M|/Af 2 . The first three terms of the Taylor 
series of II (q 2 ) are evaluated analytically up to order a 2 . The first two are required 
to derive the approximation, the third can be used to demonstrate the smallness of 
the neglected terms. The paper improves earlier results based on the leading term 
oc GfM 2 c? s . Results for the subleading contributions to Ar and the effective mixing 
angle sin 2 6 are presented. 
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1 Introduction 



The precision achieved in recent measurements of electroweak observables has sur- 
passed by far earlier expectations. The predictions for these quantities and the relations 
between them, which are based on the Standard Model (SM), are strongly affected by 
radiative corrections. Particularly important are those induced by virtual contributions 
from the heavy top quark 0. Their verification provides an important test of the the- 
ory and its quantum corrections. The agreement between the value of M t suggested by 
the CDF and DO collaborations § of M t = 176 ± 8 ± 10 GeV and M t = 199l£? ± 22 
GeV, respectively, and deduced indirectly from precision measurements of M t = ^^-H-lo 
GeV [|IJ, constitutes a triumph of the SM and a verification of the quantum corrections 
with increasing precision. A more refined understanding of these effects is on the agenda, 
including two- and even the dominant three-loop contributions [[|, || . 

As far as contributions from the top-bottom multiplet are concerned, the perturbative 
results for the self-energies Uww{q 2 ), ^-zz(q 2 ), n 77 (g 2 ) and Il 7 ^(g 2 ) are available in the 
literature for arbitrary top and bottom masses up to order a s || 0, || . 

This allows evaluation not only of the leading corrections, which are governed by the 
p parameter and increase with M 2 , but also of the subleading terms. These are required 
for a complete calculation of Ar (entering the relation between Gp, M^, M§ and a) or 
of the effective mixing angle sin 2 (governing asymmetries in Z production and decay). 
Recently also three loop QCD corrections to the p parameter of 0(GpM 2 a 2 ) have been 
calculated 0, |5|. These, in turn, control the dominant terms of order GpM 2 a 2 in Ar and 
sin 2 9. 

The technique described in j|, |5| can be employed to obtain also the Taylor series 
coefficients of II(g 2 ) around q 2 = 0, in principle to arbitrary orders in q 2 and in second 
order in a s . It will be demonstrated below that the two lowest terms in the expansion of 
Ar in M§ /M 2 provide an excellent approximation to the full answer in one- and two-loop 
approximation, corresponding to a° and a] corrections. This justifies the expectation 
that also in order a 2 the leading terms oc GpM 2 a 2 plus the subleading terms provide 
an adequate description of the complete answer for Ar and sin 2 G. This is verified by 
calculating the terms oc M|/M 2 which indeed turn out to be negligible. The complete 
result for Ar and sin 2 G to order a 2 is therefore at hand. 

2 The M\y — Mz connection and the effective mixing 
angle 

It has become customary to express the magnitude of radiative corrections in the relation 
between M w , M z , Gp and a through the quantity Ar, defined through 
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The influence of a heavy quark doublet can be expressed through the transversal parts of 
the gauge boson self energies: 
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where Il 77 (g 2 ) = Il 77 (q ,2 )/g 2 . In terms of the transversal parts of vector and axial current 
correlators H v and U A the building blocks for Ar are given by 
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with Vi = 2/3 — 4s 2 Qi. The sin and cos of the weak mixing angle are denoted by s and 
c. The "singlet" contribution to the axial current correlator U A ' S which originates from 
double triangle diagrams occurs first in order a 2 and has been displayed separately. 

Throughout this paper the mass of the bottom quark will be neglected. To circumvent 
the mass singularity which originates from the bottom loop contribution I1 77 (0) is replaced 
by Il 77 (g 2 = M§ ). The difference is accounted for by dispersion relations with input from 



the actual measurement of a(e + e~ — > hadrons) in the low energy region [ 1Q[ . 

The evaluation of Il(g 2 ) up to three-loop diagrams involving massless diagrams only 
is performed with the FORM [|ll|] package MINCER [12| implementing an algorithm 
developed in |]13| . For diagrams involving a heavy top quark the approximation based on 
a Taylor series of Il(g 2 ) around q 2 = 0, which is equivalent to an expansion in M\jM x 
leads to an adequate approximation: 



t ■ 



Af tb = Ar tb - fl 77 (0) + Refl 77 (Mf) 
c 2 /IW(0) U zz (0)\ 



+- 



+ 



M? 




w 



WW 



Ml 



(o) - n: 



zz 



Re 



^zz{M\% 



Ml 



+ n 



Ml 



77 



+ Refl 77 (M|)| f) -n / w (0) 



n' (o)| t Ml-n'^(o) 



M\ 



w 



(4) 



The first line in this expansion is determined by the p parameter, leads to terms oc GfM 2 , 
and has been recently calculated up to order a 2 . 
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The second line leads to terms of order GpM\ and GpMf In M 2 jM\ and will be 
calculated below up to order a 2 . The subsequent terms are of order GpM|/M 2 . Even 
after charge and mass renormalization H(q 2 ) exhibits (in dimensional regularization) 1/e 
singularities which cancel in the combination ([|). 

The reliability of the M 2 Z /M 2 expansion is demonstrated in Fig. [I] for the one- and two- 
loop results respectively. The full answer (solid line) is compared to the approximation 
based on the term quadratic im M t (dotted line) and the approximation including constant 
plus In M z /M 2 terms (dashed line). The latter provides an excellent approximation in the 
range M t > 150 GeV; corrections of order GpMf/M 2 may safely be ignored (dash-dotted 
line) . 

The leading terms of the expansion are given by WM (X t = GfM 2 / (8\/2tt 2 )) 
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The full result is given e.g. in |7|, [L4|. 

A second quantity of practical interest is the effective weak mixing angle, which governs 
in particular the asymmetries [l| in Z boson production and decay. It is related to 
sin 2 = 1 — .\/f r /.\ l z through a correction factor 



sin 2 6 = (1 + Ak) sin 2 9 
which in turn is influenced by the polarization functions of Eq. 
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The full analytic result and the approximation based on the Mz/Mf expansion are com- 
pared in Fig. |2], adopting the same notation as in Fig. [TJ. From these figures and from 
Eqs. (H,^) it is evident that (for M t « 180 GeV) the next-to-leading corrections amount 
to about 25% of the GpM 2 terms. The next-to-next-to- leading terms of order GpM\jM 2 
are below 1.2%. These considerations justify the restriction of a 2 corrections to the first 
two or at most three terms in the Mf/M 2 expansion. 



3 Three-Loop Corrections 

In addition to the massless diagrams there are two further types of integrals: The evalu- 
ation of the derivatives of IT(g 2 ) resulting from diagrams where the massive top quark is 
coupled to the W or Z is reduced to the evaluation of tadpole integrals discussed in |Q, ||] . 
The derivatives are obtained through Taylor expansion of the respective integrands up to 
0(q 2 ) 

n^(g) = ^JI(g 2 ) + q,q u Tl 2 (q 2 ) (8) 

and projecting out the transverse part. Following [|T3| , |16f the resulting tadpole integrals 
are subsequently reduced to the master set listed in [f|, ||. Diagrams involving external 
massless quark loops and internal top loops (and the anomaly graph) are treated with 
the large mass expansion technique JT7 . 



Setting Ca = 3, Cp = 4/3 and /i = m 2 a fairly compact form for the subleading parts 
of Af t b and An t b is obtained (x t = Gpfh 2 / (8 s/2-k 2 ), T z = lnMf/m 2 ): 
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The formula for 5#j^g can be found in 



This result is formulated in terms of the MS coupling a s and the mass m t . Employing 
the two-loop relation between the MS mass and the pole mass M t the MS results are 
easily expressed in terms of M t . The final result, after setting rif = 6, reads (/i 2 = M 2 ): 



Af 



OS 
tb 



-^2,X t 1 + 



Ml 



(0.3333 - 0.6667 l z + (-0.5926 + 0.8889 l z )s' 



+ 



- 0.4 + 1.144s 2 - s 4 



5 



+- 



a x r 



7T 



2.8599 + tt|( - 1-5640 - 0.6667/ z + (0.1022 + 0.8889 l z )s' 



+ 



7T 



'M 



2\ 2 



ill 



1.312 + 3.573s 2 - 3.582s 4 



" " - 14.594 + — f f - 17.224 + 0.08829 / z + 0.4722 Z 



+ (22.6367+ 1.2527 l z - 0.8519 / z 2 )s 2 + 



'Mf 



- 7.7781 - 0.07226 Z z 



+ 0.004938 1£ + (21.497 + 0.05794 l z - 0.006584 l 2 )s 2 - 21.0799s 



OS 
tb 



-%3X t 



1 + 



( - 0.3333 + 0.6667/ z + (0.2963 - 0.4444/ z )s 2 ) 



+ 



'Ml 



(0.4 - 0.8222s 2 + 0.5s 4 ) 



+ - 



M~ 

2.8599 + — f (1.564 + 0.6667 l z + (-0.051103 - 0.4444 l z )s 2 ) 



M' t 



'Ml 
Ml 



(1.312 - 2.682s 2 + 1.791s' 



+ 



Mi 



14.594 + —|f 17.224 - 0.08829 l z - 0.47222/ 

Mi v 



+ (-11.3184 - 0.6263 l z + 0.4259 Zj)s 2 ) + 



M 2 



7.7781 + 0.072263 Z z 



- 0.004938 1% + (-16.0186 - 0.02897 l z + 0.003292 /j)s 2 + 10.54s 4 ) 



+i 



■ Ml 
M? 



- 1.396s 2 + 1.862s 4 + —(-1.396s 2 + 1.862s 4 

7T 



( — 

7T 



+ 



M 2 



1.968 + 2.676 l z )s 2 + (2.6235 - 3.5682 l z )s 4 
2 ((-0.09102 + 0.02069 l z )s 2 + (0.1214 - 0.02758 l z )s 



where l z = In Ml /M 2 . In these formulae also the terms of order GFOt 2 M z /M 2 are given. 
As shown in Fig. [I] and and Table [I] their effect of the numerical result is extremely 
small and can safely be neglected. 

In the lowest diagrams of Fig. [I] and [2] the 0(a 2 ) corrections of Af#, and A% are 
presented as functions of M t and the quality of the M|/M 2 expansion is confirmed. The 
difference between the quadratic term (dotted line) and the constant plus log term (dashed 
line) amounts to about 25%. Adding the subsequent term proportional M|/M 2 (dashed 
dotted line) barely affects the answer. The numerical effects on My/ and sin 2 are given 
in Table |1|. The numbers are obtained with the following input data: a s (M 2 ) = 0.1092 
(corresponding af\M 2 z ) = 0.120), M t = 175 GeV, M z = 91.188 GeV, M H = 300 GeV, 
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M 2 


Mf + const. 


M 2 - 


h const. + 1/M 2 


5M W /M W (OS) 


0.00677 


0.00838 




0.00825 


5 sin 2 ©/ sin 2 © (OS) 


-0.01618 


-0.01832 




-0.01814 


5M W /M W (MS) 


0.00674 


0.00833 




0.00821 


5 sin 2 ©/ sin 2 © (MS) 


-0.01610 


-0.01823 




-0.01804 



Table 1: Numerical results for the M 2 , the constant plus logarithmic and the 1/M 2 con- 
tributions. For the numerical evaluation of 5 sin 2 ©/ sin 2 © only the real part of Ak^ is 
taken. 



5M W in MeV 


a 


ct 1 


ct 2 


M 2 


611.9 


-61.3 


-10.9 


const. 


136.6 


-6.0 


-2.6 


1/M 2 


-9.0 


-1.0 


-0.2 



Table 2: The change in Mw separated according to powers in a s and M t in the on-shell 
scheme. 



ct = 1/137.04 and Gp = 1.16639 10~ 5 GeV~ 2 . In each column the terms of order ct°, a] 
and ct 2 belonging to the corresponding expansion in the top mass are added up. For Ar 
in Eq. (H) we used Ar = Act + Af tb + S rem with Act = 0.05940 Hl0| . 5 rem contains all 
contributions of order ct which are not present in the other two pieces and can e.g. be 
found in fl9 |. In Table |^ the contributions are listed separately according to powers of a s 
and M t . One observes that the absolute prediction for the W mass is changed by —10.9 
MeV if the ct 2 M 2 term is added to the full two-loop result. This increases to —13.7 MeV 
if also the constant and the 1/M 2 suppressed terms terms are added. (For a fictitious top 
mass of 100 GeV the numbers would be —4.2 MeV and —7.9 MeV respectively.) 

Summary: Top mass dependent corrections to relations between electroweak parame- 
ters have been calculated up to order ct 2 , with the help of an expansion in M|/M 2 . The 
quality of the approximation has been confirmed in one-, two- and three- loop approxima- 
tions. 
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Figure Captions 



Figure 1: Af t & as a function of M t . The dotted, dashed and dash-dotted curves correspond 
to an increasing number of terms in the approximation. In the figures a value s 2 = 0.2321 
is chosen. 

Figure 2: AKtb as a function of M t . The conventions are the same as for Af t b- Here only 
the real part is plotted. 
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